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We consider the Casimir interaction between a ferromagnetic and a non-magnetic mirror, and show 
how the Casimir effect gives rise to a magnetic anisotropy in the ferromagnetic layer. The anisotropy 
is out-of-plane if the non-magnetic plate is optically isotropic. If the non-magnetic plate shows a 
uniaxial optical anisotropy (with optical axis in the plate plane), we find an in-plane magnetic 
anisotropy. In both cases, the energetically most favorable magnetization orientation is given by the 
competition between polar, longitudinal and transverse contributions to the magneto-optical Kerr 
effect, and will therefore depend on the interplate distance. Numerical results will be presented for 
a magnetic plate made out of iron, and non-magnetic plates of gold (optically isotropic), quartz, 
calcite and barium titanate (all uniaxially birefringent). 



I. INTRODUCTION 

Even 60 years after its discovery, the Casimir effect 
continues to fascinate [l|. Originally calculated as an 
attractive force between two neutral metallic plates due 
to the zero-point energy of the electromagnetic field, the 
Casimir effect nowadays plays a role in fields as diverse 
as gravitation, cosmology, atomic and molecular physics, 
and quantum field theory [2|]. From a condensed matter 
perspective, research has been mainly preoccupied with 
the inclusion of several contributions that are present in a 
real experimental setup, e.g., boundary roughness, shape 
and geometry of the mirrors, non-zero temperature, and 
dielectric properties of the boundaries (a thorough review 
can be found in Ref. 0). With the theory in place, sev- 
eral recent experiments have verified the existence of the 
Casimir force @, 1, H, H, 0, B ©, ) . 

Whereas most publications deal with optically 
isotropic mirrors, the Casimir effect between anisotropic 
mirrors is a rather recent topic. The van der Waals in- 
teraction between uniaxial anisotropic media of infinite 
thickness has been considered in Ref. [ll[ and formed the 
starting point for predictions of a torque between two 
uniaxial birefringent plates [HI, 0, [3] • The zero-point 
energy depends on the relative orientation of the plates 
which gives rise to a torque that seeks to align their prin- 
cipal axes. Similar effects are present for plates in which 
the optical anisotropy is induced by corrugations [l5|. 
Furthermore, in a search for repulsive instead of attrac- 
tive Casimir forces, and stimulated by the advent of so- 
called magneto-dielectric metamaterials (which have a 
non-negligible magnetic response on top of the electric 
one), the effects of anisotropy of the magnetic permeabil- 



ity on the Casimir effect were studied in Refs. [16j, Il7j . 
The Casimir force between ferromagnetic mirrors, which 
of course cannot be described by a magnetic permeability, 
has been treated by us in a recent publication [l8|, [l9| . In 
this case, the boundary conditions for the zero-point elec- 
tromagnetic field are influenced by the magneto-optical 
Kerr effect, and we have found that this gives rise to a 
magnetic contribution to the Casimir force. 

In the present paper, the force between a single ferro- 
magnetic mirror opposite of a non-magnetic one is con- 
sidered. Interestingly, we find that the Casimir effect will 
give rise to a magnetic anisotropy in the ferromagnet. We 
will first study a setup where the non-magnetic mirror is 
optically isotropic, and observe an out-of-plane magnetic 
anisotropy in the ferromagnetic layer. For this setup, the 
anisotropy energy arising from the Casimir effect is or- 
ders of magnitude smaller than the surface anisotropy 
and magnetocrystalline anisotropy. Therefore we will 
propose another scheme where the non-magnetic mirror 
is made of a uniaxially birefringent material (with optical 
axis in the plane). The Casimir effect then induces an 
in-plane magnetic anisotropy that can be detected more 
easily. In both cases, the magnetic orientation with the 
lowest energy results from a competition between several 
contributions to the magneto-optical Kerr effect, and will 
depend on the distance between the plates. 



II. GENERAL THEORY 

The Casimir energy (per unit area) between two plates 
A and B can be expressed in terms of the reflection co- 
efficients of the plates as 0, [2(| 
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where Tilo is the electromagnetic field frequency, k± is the 
component of the wavevector perpendicular to the mir- 
rors, and if is the angle between the incidence plane and 
the X-axis. Furthermore, R 



A/B 



are 2x2 matrices con- 



taining the reflection coefficients of the plates (evaluated 
at imaginary wavevector and frequency) : 
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The index s (resp. p) corresponds to a polarization with 
the electric field perpendicular (resp. parallel) to the 
incidence plane. We adopt the convention that the p 
axis remains unchanged upon reflection. 

In this paper, we will consider one of the plates to be 
ferromagnetic. The reflection coefficients are then influ- 
enced by the magneto-optical Kerr effect, and depend on 
the magneto-optical constant Q. Since Q is small, we will 
only keep terms up to first order in Q, in which case the 
reflection matrix can be written as R = Rq + AR, with 



of the diagonal (e xx ) and off-diagonal (e xy ) components 
of the dielectric tensor [13, \2M ■ 
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where we have defined: £ = \J oj 2 (e xx (iu)) — 1) + (fcj_c) 2 . 



III. OUT-OF-PLANE MAGNETIC 
ANISOTROPY 
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depends on the magnetization. The magnetization ori- 
entation has been written out explicitly in terms of the 
spherical angles 9 and <f> (the Z-axis points into the fer- 
romagnetic plate) . The terms r ss and r pp are the Fresnel 
coefficients for a normal metallic mirror, while the terms 

r sp> r ^p> an d Ar pp arise from the polar, longitudinal and 
transverse Kerr effect respectively and are due to the 
magnetization. All coefficients can be expressed in terms 



Using the reflection coefficients above, we will now con- 
sider the Casimir effect between a non-magnetic isotropic 
mirror A and a ferromagnetic mirror B. For this partic- 
ular case, mirror A is described by a reflection matrix 
= R^i an d mirror B by Rb = + AR as defined in 
Eqs. ©, ([3]) and |5]). Making an expansion in A R of the 
logarithm in Eq. (JT|) and keeping only terms in the lowest 
non- vanishing order, one obtains the following expression 
for the Casimir energy: 
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In this expression, £q is independent on the magnetiza- 
tion orientation and describes the non-magnetic contri- 
bution to the Casimir energy, while 
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are terms arising respectively from the polar, longitudinal 
and transverse contributions to the magneto-optical Kerr 
effect. 

The Casimir energy in Eq. ((SJ) is of course indepen- 
dent of the in-plane magnetization, because the consid- 
ered system is invariant with respect to rotations around 
the normal to the plates. However, the energy depends 



on the angle 9 between the magnetization and the nor- 
mal to the plates, i.e., the Casimir effect gives rise to an 
out-of-plane magnetic anisotropy. It is important to note 
that the contribution coming from the polar Kerr effect, 
£ is proportional to cos 2 9 and thus tends to align the 
magnetization parallel to the plates (9 = ir/2 or 3ir/2), 
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while the longitudinal and transverse contributions £{ 2 
are proportional to sin 2 9 and have their energy mini- 
mum when the magnetization is perpendicular to the 
plates (6 = 0, it). The physical parameters (dielectric 
constants) of the plates, and as we will see, the interplate 
distance, will decide which of these competing terms is 
dominant. 

We have numerically calculated the different contribu- 
tions for a setup where the non-magnetic plate is made of 
gold, and the magnetic plate is made of iron. The reflec- 
tion coefficients in Eq. ([5|) will depend on the dielectric 
tensor of these metals. In order to take the interband 
transitions properly into account, we have calculated the 
dielectric tensor from the available optical data, instead 
of reverting to a simple plasma or Drude model. The di- 
agonal part of the dielectric response function e xx (uj) of 
Au and Fe, evaluated at real frequencies, is tabulated in 
Ref. [12] • A significant amount of data points is available 
in the range 1(T 1 . . . 10 4 eV for Au and ft)" 3 . . . 10 4 eV for 
Fe. At lower frequencies, where interband transitions do 
not play a significant role anymore, we have extrapolated 
the data with a Drude model using a plasma frequency 
hw p = 9 eV and a relaxation time H/t = 35 ■ 1 -3 eV 
for Au, while faj p = 3.54 eV and H/t = 19 • ICT 3 eV 
for Fe [23|. Optical data for the off-diagonal element of 
the dielectric tensor e xy (uj) of Fe is rather scarce. We 
have used data points between 0.1 eV and 6 eV found in 
Ref. (23|, and assumed e xy {oj) = out of this range. The 
dielectric constants e xx (\oj) and e xy (ioj) at imaginary fre- 
quencies, which are needed for evaluating the reflection 
coefficients in Eq. ([5]), were then obtained by using the 
Kramers-Kronig relations. Details of the whole proce- 
dure, together with plots of the resulting dielectric ten- 
sors can be found in Ref. [13] and our earlier work [Tgj | . 

Finally, the Casimir energy per unit area is calculated 
by numerical integration of Eqs. J7J), for interplate dis- 
tances D ranging from I nm to 5 fim. Results for the 
Casimir energy amplitudes £^ , £ {' , and £ J are given in 
Fig. Ha). The corresponding force contributions are ob- 
tained by deriving the energy expressions with respect to 
the interplate distance D: 

and are shown in Fig. QJb). 

For distances larger than approximately 10 nm, corre- 
sponding to the experimentally relevant regime, the con- 
tribution £ arising from the polar Kerr effect is clearly 
dominant. Looking back at Eq. ([5]), this means that 
the lowest energy state at these distances is one where 
the magnetization is perpendicular to the plates. How- 
ever, at smaller distances there is a crossover to a regime 
where the transverse contribution &\ dominates, and con- 
sequently the preferential magnetization direction will in- 
terestingly change from perpendicular towards parallel to 
the plates upon decreasing the interplate distance. Look- 
ing at the distance dependence of the force (per unit area) 
in Fig. [TJb), we have found that for long distances the 
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FIG. 1. Polar, longitudinal, and transverse contributions to 
the Casimir energy per unit area (a) and the Casimir force 
per unit area (b) between an Fe and an Au mirror. 

polar contribution varies like D~ 13 / 2 , compared to D~ 8 
for the other contributions. For very short distances, 
the transverse force varies like D~ 3 , while the polar and 
longitudinal terms show (approximately) a depen- 
dence. 

In order to estimate the Casimir force in a typical ex- 
periment, we consider two disks of radius R — 10 /im, at a 
distance D = 100 nm from each other. For this setup, we 
obtain a force amplitude F = J-ttR 2 on the order of a few 
fN. The anisotropy could thus possibly be detected with 
state-of-the-art mechanical (force) measurement schemes 
(aN resolution is reached in magnetic resonance force 
microscopy [III H(| Ht||), but presumably not via some 
magnetic measurement because the anisotropy contribu- 
tion due to the Casimir effect is considerably smaller (by 
several orders of magnitude) than the shape anisotropy 
and the magnetocrystalline surface anisotropy of the Fe 
plate. We have therefore looked at a different configu- 
ration where the latter anisotropies are essentially zero. 
In particular, we considered a setup where the magnetic 
mirror is now placed opposite of a non-magnetic mirror 
with uniaxial (optical) anisotropy. The latter can, e.g., 
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be made from a uniaxial birefringent material with the 
optical axis in the plane. For this kind of system, we have 
found that the Casimir effect will create an in-plane mag- 
netic anisotropy, which could be much easier to detect. 
Our results will be presented in the next section. 



IV. IN-PLANE ANISOTROPY 

The reflection coefficients of a mirror with uniaxial op- 
tical anisotropy can be expressed in terms of the ordi- 
nary (e ) and extraordinary (e e ) dielectric constants of 
the anisotropic material. The corresponding equations 
can be found, e.g., in Refs. [H,[29j]. We will consider a 
mirror with the optical axis in the plane, for which one 
obtains (at imaginary frequency and wavevector) [29j : 
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Cq = \J (e - a 2 ) + (afc^c) 2 , 

and where a = sin£ and (3 = — cosC, with £ the angle 
between the optical axis and the X-axis. The dielectric 
functions e and e e are evaluated at imaginary frequen- 
cies. Please note that in comparison with the expressions 
in Ref. [29| , both r pp and r ps have a minus sign. This re- 
sults from a different convention in Ref. (29|, where the 
p axis reverses si gn u pon reflection, whereas we keep this 
axis unchanged [301 ] . Furthermore, the definition of a 
and P is different because we assume the incidence plane 
to be the OYZ plane as in Ref. (20| . whereas this is the 
OXZ plane in Ref. [H- 

It is clear that the angular dependence of these reflec- 
tion coefficients is much more complicated than what we 
had for a magnetic mirror in Eq. ([4]) . Therefore it is more 
difficult to make analytical progress, and we have opted 
for a direct numerical evaluation of the Casimir energy in 




co (eV) 

FIG. 2. Dielectric functions for quartz (a), calcite (b), and 
barium titanate (c). Solid and dashed lines correspond to the 
ordinary and extraordinary dielectric function respectively. 



Eq. fl}. We have considered a setup where one ferromag- 
netic (Fe) plate with magnetization in the mirror plane is 
opposed to a non-magnetic uniaxially anisotropic mirror. 

The dielectric tensor of the Fe plate is calculated di- 
rectly from optical data, as has been described in the 
previous section. The uniaxial plate is assumed to be 
made of one of the birefringent materials quartz (Si02), 
calcite (CaCOs) or barium titanate (BaTiOs). Since op- 
tical data for these materials is not available in a large 
enough frequency range, we have used a two-oscillator 
model to describe their dielectric function. Although this 
model is quite simple, it gives a rather good description 
of the dielectric tensor in inorganic materials [30L l3lj , and 
has been used already to estimate the Casimir torque be- 
tween two birefringent plates [l3:] . Within this approach, 
both the ordinary and extraordinary dielectric function 
of the uniaxial mirror will be modeled by: 



e(iw) = 1 + 
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(11) 



where uir(uv) an d ^IR(UV) are the characteristic ab- 
sorption frequencies and absorption strengths in the in- 
frared (ultraviolet) range. These constants were calcu- 
lated in Refs. [3l|,[33] from the available optical data, and 
were summarized for the ordinary and extraordinary di- 
rections of quartz, calcite and barium titanate in Table 
I of Ref. [IH . The resulting dielectric functions are plot- 
ted in Fig. [2] [33[. From this figure, one sees that quartz 
shows the weakest birefringence and has e e > e for all 
frequencies. Barium titanate is the strongest birefringent 
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FIG. 3. Amplitude of the anisotropy energy (per unit area) 
for a magnetic Fe mirror placed opposite of a non-magnetic 
birefringent mirror made out of quartz, calcite or barium ti- 
tanate. Interplate distance D — 10 nm (a) and D = 5 fim (b). 
The amplitude is multiplied by 10 for quartz in both figures, 
and for calcite in (b). 

material and has e e < e Q , whereas calcite is somewhere in 
the middle and shows e e > e for small frequencies and 
e e < e for larger frequencies. 

We have calculated the Casimir energy £ (A(f) as a 
function of the angle A(f> between the magnetization of 
the magnetic mirror (assumed to be in-plane) and the 
optical axis of the non-magnetic one. Results are shown 
in Fig. [3] for the three considered birefringent materials, 
and at distances of 10 nm and 5 /im. We have plotted the 
quantity AS = S (A</>) — S m i n so that the energy minimum 
will always lie at zero. All curves could be fitted very 
well with either a sin 2 A(f> or a cos 2 A(f> dependence, as 
depicted by the solid line connecting the calculated data 
points. 

We observe that at the shortest distance (10 nm in 
Fig. GDJa)), the curves for BaTiC>3 and CaCC>3 show a 
cos 2 A(j) behavior and thus favor a magnetization that is 
perpendicular to the optical axis, whereas the curve for 
SiC>2 varies like sin 2 Acf) and tends to align the magne- 
tization with the optical axis. This difference is due to 
the different sign of the quantity Ae = e e — e G for the 
three materials in Fig. [2l and we find that the energy 
minimum corresponds to a situation where the magneti- 
zation is oriented perpendicular to the dielectric axis with 
the largest value of e. One should note that although Ae 
changes sign at some intermediate frequency for calcite, 
the higher frequencies (where Ae > 0) will dominate the 
integrals in Eq. (TI| at small distances, and calcite there- 
fore shows the same angular dependence as BaTiC>3 in 
Fig. Ha). 

Interestingly, at a longer distance (5 ^m in Fig. (3^b)), 
the angular dependence is exactly opposite, at least for 
quartz and barium titanate: we have a sin 2 A</> behavior 
for BaTiC>3 and a cos 2 A</> for quartz (see Fig. [3j so that 
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FIG. 4. Amplitude of the anisotropy energy (per unit area) as 
a function of interplate distance D for a magnetic Fe mirror 
placed opposite of a non-magnetic mirror made out of quartz, 
calcite, or barium titanate. 

now the energy is minimized for a magnetization parallel 
to the dielectric axis with the largest value of e. In the 
previous section, we observed a similar change of the an- 
gular dependence for a setup with an optically isotropic 
mirror, and explained it by a competition between the 
polar and transverse Kerr contributions. Although the 
polar effect is zero in the setup we consider now (the 
magnetization lies in the plane, so 9 = tt/2 in Eq. Q), 
the two remaining longitudinal and transverse contribu- 
tions can still compete with each other and will lead to 
the observed angular crossover. It is however not possi- 
ble to address both contributions separately, since mixed 
terms of the form rl p Ar pp will also appear in the energy 
expression. 

Calcite is a special case, as it has the same cos 2 A(j> 
behavior both in the limit of short and long distances 
(see Fig. [3]). In principle, the competition between lon- 
gitudinal and transverse Kerr effect contributions would 
lead to the same crossover as in the other two considered 
materials. However, at long distances the low frequen- 
cies will determine the integrals in Eq. (p}. As depicted 
in Fig. HJb), one has Ae < at these low frequencies, 
compared to Ae > for the higher frequencies which de- 
termined the integrals at short distances. It is this sign 
change of Ae that will cancel the angular dependence 
crossover due to the Kerr effect competition. Neverthe- 
less, there will be an intermediate distance regime where 
the higher frequencies with Ae > are still the most im- 
portant, and the competition between the Kerr contri- 
butions is decisive in changing the angular dependence 
from cos 2 Atp to sin 2 A(f>. 

This becomes clearer in Fig. [4j where the amplitude of 
the A<fi oscillations is shown. Both the curves for BaTiC>3 
and quartz show a single kink. At the kink, the sin 2 A</> 
and cos 2 A</> contributions to the anisotropy energy are 
of equal magnitude so the net energy amplitude goes to 



6 



zero. This kink thus signals a crossover point, and as ex- 
plained above, one clearly observes two such crossovers 
for calcite but only a single one for quartz and barium ti- 
tanate. Furthermore we see from Fig.[3]that, as expected, 
the anisotropy energy is proportional to the amount of 
birefringence of the considered material, and thus largest 
for BaTiC>3 and smallest for quartz. 

An experimental detection of these anisotropy effects 
will be quite challenging since the calculated energies are 
rather small. However, the studied system has the ad- 
vantage to lend itself to different measurement setups. 
Firstly, a force measurement could be done by evaporat- 
ing a droplet of iron on a cantilever and placing it over 
a birefringent plate. A resonant rf-field can be used to 
rotate the magnetization, and the Casimir force will re- 
sult in oscillations at the resonance frequency of the can- 
tilever. For an iron droplet with a curvature radius R = 
100 /im, the force amplitude F can be estimated by using 
the force proximity theorem: F — 2irRA£. At a distance 
of lOnm, the force is F = 60 aN for BaTiC>3. Force detec- 
tion with aN-resolution has already been shown in similar 
magnetic resonance force microscopy setups [25l . l26l |27| . 
Secondly, the anisotropy energy will give rise to a me- 
chanical torque dAE/dAcj) that seeks to rotate the plates 
towards their minimal energy position. For plates with 
a radius of R = 100 /im, the torque at an interplate 
distance D — 10 nm is « 10~ 21 Nm. A possible s etup 
to demonstrate such torques can be found in Ref. [13 |. 
where both plates are immersed in a liquid one on top 
of the other. With a suitable selection of material pa- 
rameters, a repulsive Casimir force will balance gravity 
and a birefringent disk can be held floating above a hard 
ferromagnetic plate. The Casimir torque will then ro- 
tate the plates back to their minimum energy orientation 
once they are brought out of their equilibrium position. A 



third measurement scheme would consist of detecting the 
magnetic anisotropy energy of the ferromagnetic plate di- 
rectly with a magnetic measurement. By using a cylindri- 
cally shaped ferromagnetic mirror cut along the correct 
axis, both the shape and magnetocrystalline anisotropy 
can be made small. Furthermore, no anisotropy contri- 
bution intrinsic to the ferromagnetic plate will depend on 
the orientation of the birefringent plate, so that parisitic 
effects can in general be filtered out efficiently. Finally, 
we would like to mention that instead of using a bire- 
fringent material such as barium titanate to create the 
uniaxially anisotropic mirror, the optical anisotropy can 
also be realized by means of some appropriate nanostruc- 
turing, e.g., by creating parallel stripes on an isotropic 
gold mirror. The optical anisotropy can then be made 
much larger, thereby increasing the magnetic anisotropy 
effects (a factor of 1000 is possible [Hj]). A detailed cal- 
culation for this case is not at all straightforward, and 
beyond the scope of the present paper. 

V. CONCLUSION 

We have considered the Casimir effect between a fer- 
romagnetic and a non-magnetic mirror. If the non- 
magnetic mirror is optically isotropic, the Casimir ef- 
fect is shown to give rise to an out-of-plane magnetic 
anisotropy in the ferromagnetic plate. For a non- 
magnetic mirror with uniaxial anisotropy, we found an 
in-plane anisotropy. In both cases, the angular position 
with the lowest energy depends on the distance between 
the plates and is determined by the competition between 
different contributions to the magneto-optical Kerr ef- 
fect. Different setups for measuring the effect have been 
proposed. 
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